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Abstract 

We review the use of superconductors as a playground for the experimental study of front rough- 
ening and avalanches. Using the magneto-optical technique, the spatial distribution of the vortex 
density in the sample is monitored as a function of time. The roughness and growth exponents 
corresponding to the vortex 'landscape' are determined and compared to the exponents that char- 
acterize the avalanches in the framework of Self-Organized Criticality. For those situations where a 
thermo-magnetic instability arises, an analytical non-linear and non-local model is discussed, which 
is found to be consistent to great detail with the experimental results. On anisotropic substrates, 
the anisotropy regularizes the avalanches. 
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FIG. 1: Vortex density (along the vertical axis) in a thin YBa2Cu30y- superconductor as a function 
of position. Note the similarity with a pile of sand, including the surface roughness. 



I. INTRODUCTION 

In this review we discuss the penetration of magnetic flux in thin type-II superconduct- 
ing slabs subject to an increasing magnetic field oriented normal to the slab. In type-II 
superconductors, magnetic flux penetrates in the form of vortices, which behave as repulsive 
particles moving in the plane of the slab. The vortices can only enter (or leave) at the edge 
of the slab. Their motion, however, is hindered by pinning (due to disorder in the sample), 
which leads to a gradient in vortex density, somewhat similar to the slope of a sand pile, 



see Fig. [TJ This analogy was already noted a long time ago by de Gennesjlj]. In the present 
paper, however, this analogy is extended to the surface roughness of the pile and to the 
avalanches occurring on its surface. 

The growth and roughness exponents of this rough surface are determined. Also the 
properties of the avalanches occurring on the surface are measured and it will be shown that 
they have a power law size distribution function and obey "finite size scaling". In addition, 
we find that exponent scaling relations of Paczuski, Maslov and BakQ] are obeyed to reason- 
able accuracy. All these observations together make a strong case for the vortex pile system 
to belong to the class of Self-Organized Criticality (SOC) 3] systems. To investigate the role 
of disorder, we also study thin films of Niobium, which are slowly subjected to a hydrogen 
atmosphere. Hydrogen is absorbed by Niobium and locally destroys the superconductivity, 
thus increasing the amount of disorder. We find that a minimum amount of disorder is nec- 
essary for SOC to occur. In other superconductors, the properties of huge thermo-magnetic 
avalanches are studied. These are described very nicely by a recently developed model ^| 
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and can be simulated to high degree of correspondence. 

II. EXPERIMENTAL TECHNIQUE 

Images of the 2-dimensional (i.e. in the plane of the thin film sample) density of vortices 
are made using a magneto-optical technique [sj]. On top of the sample and in immediate 
contact with it, a mirror layer and (directly above this) a 4 /im thick Yttrium Iron Garnet 
film with large Faraday effect are placed. The sample assembly is placed on a custom built 
insert with X, Y, Z, 6, </>, //-controls for proper alignment, which fits in a commercial 1 T or 7 T 
Oxford instruments cryostat. The insert is equipped with lenses, a semi-transparent mirror 
and polarizers, making it effectively a polarization microscope. Where ever there is a local 
magnetic field in the sample, the incoming polarization vector is turned, thereby increasing 
the detected intensity. Thus the intensity in the camera yields directly the magnetic fiek 
and hence vortex density for each position in the sample. Using a modulation scheme jf], E| 
this intensity map becomes linear in field and much more sensitive. 

III. ROUGHENING OF THE VORTEX LANDSCAPE 

As mentioned above, the vortex density map of the sample is generally not smooth (see 
Fig. but displays some roughness. For example, one may define a flux front as the 
'foot' of the vortex 'mountain' as indicated by the white line in Fig. ^ The exact position 
is determined by the cross section between the vortex 'mountain' and a horizontal plane 
at a level of 3 times the noise in the vortex-free region. Generally 0, for such fronts, the 
r.m.s. width w growths initially with time asw~ t 13 . After some time t x , however, lateral 
correlations extend over the whole area of the pile and w does not increase anymore: it is 
limited to a pile-size dependent value w ~ L a where L is the linear size of the pile and a is 
the roughness exponent. At t x , both relations hold from which 
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To find a and f3 to higher accuracy, we use the correlation function [8( 
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FIG. 2: Flux fronts in a YBa2Cii307-sample of 80 nm thick that was deposited on a NdGaO"3 
substrate for different applied external fields. The first front (bottom) is recorded at 1 mT and 
subsequent fronts are at 1 mT intervals. 



where 8h — h — h and h(x,t) is the height of the interface at position x and time t, 
where h = (h) x ^ 0L ^ and (•) denotes an average over r\. From this correlation function, one 
determines the growth exponent (3 and roughness exponent a using t 8] C (x, 0) ~ x a and 
C(0,t) ~t". 

For YBa2Cu307-samples of 80 nm thick that were deposited on NdGaO^ substrates, we 



recently found for the roughness of the 2- dimensional surface of the pile 

a = 0.8(1) (3) 
d = 0.6(1) 

The values of these exponents are consistent with the Q-EW (quenched- disorder 
Edwards-Wilkinson) model. Apparently the vortices move on a substrate with randomly 
distributed quenched pinning sites, as is indeed expected for a high-T c superconductor. Sim- 
ilar values for a and 8 for the flux fronts in YBa2Cu30y thin films were previously obtained 



by Surdeanu et al. 



12| . In gen- 



llj an d f° r a m Niobium thin films by Vlasko-Vlasov et al. 
eral, the values of the exponents for the 2-dimensional surface are lower than for the fronts. 
The reason for obtaining similar values for a and (3 for the 1-d and 2-d case is the higher 
number of defects in the sample used for 2-d case. 




FIG. 3: Size (expressed in units of the flux quantum <I>o) of vortex avalanches in a YBa2Cu307- 
sample of 80 nm thick that was deposited on a NdGaO"3 substrate as a function of time for 8 
different experiments (separated by the vertical arrows). 

IV. SMALL VORTEX AVALANCHES 

From a closer inspection of the growth process of the vortex 'mountain', shown on Fig. 
[Q it is clear that the formation of this 'mountain' and hence the flux penetration process is 
not continuous, but takes place in the form of avalanches. In our experimental procedure, 
the external field is increased in steps of 50 /iT and after each such increase an image of the 
vortex density in the sample is obtained. The size of each avalanche is calculated from 

s= l -j \AB z (v)\dv (4) 

where AB Z (r) is the difference between two consecutive images minus the 50 fiT increase 
of the external field. The resulting time series of avalanches is shown in Fig. 01 

The size distribution of the avalanches is shown for differently sized windows of observa- 
tion in Fig. E| as a log-log plot. From the straight line behavior in the limit of large systems 
(viz. large windows of observation), it is clear that the avalanche sizes are power law dis- 
tributed with P (s) ~ s~ T , just like earthquakes and many more natural phenomena^. 

The ubiquitous occurrence of power laws in many natural phenomena lead Bak to the idea 
of Self-Organized Criticality[3], a model where slowly driven systems organize themselves 
automatically to be in a metastable state where small disturbances can lead to a response of 
arbitrary size. The archetypal example is the sand pile with its critical slope (metastability) 
and the avalanches of sand sliding down its surface. Within the SOC framework (and 
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FIG. 4: Number P(s,L) of avalanches of size 5 in a YBa2Cu307-sample of 80 nm thick that was 
deposited on a NdGaC>3 substrate occuring in a window of observation of linear dimension L. The 
straigt line indicates the powerlaw in the limit for large L. 

otherwise too), the volume V of an avalanche with fractal dimension D and linear extension 
I is given by V = £ D . Since the linear size of the avalanche is limited by the system size 
L (by neccesity £ < L), one expects that the finite size effect-induced deviations from the 
straight line in Fig. i S cale a S L". Indeed the finite S1Z e scaling (FSS) plot of Fig. Eshows 
a nice data collapse and yields[9] r = 1.29 (2) and D = 1.89 (3). 



Apart from the FSS and the power law behavior there is another check one can make for 
SOC. This is based on the exponent relations of Paczuski, Maslow and Bak|2j. We discuss 
and use two such relations. 

The first scaling relation is found by calculating the volume V of the fractal avalanche. 
By definition, V ~ L D . Also, this volume equals the fractal surface times the fractal height. 
The fractal surface area is the projection of the avalanche cluster on an average (flat) plane 




FIG. 5: Finite size scaling plot of the data in Fig. |I] 
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P through the pile surface. This area is L dB where ds is the surface fractal dimension of 
the avalanche cluster. The height or thickness of an avalanche is obtained by subtracting 
the heights (in the direction perpendicular to P) of the piles before and after the avalanche. 
Since the avalanche modifies the pile only locally, this fractal thickness scales as the surface 
roughness and is proportional to L a . Hence V ~ L dB L a . Combining both expressions for V 
yields L D ~ L dB L a from which we obtain the scaling relation 

D = d B + a (5) 

The second scaling relation follows from the fact that the deviation from power law 
behavior seen in Fig. |U occurs because avalanches above a certain size 'feel' that the pile is 
finite. 
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FIG. 6: The scaling relation D (2 — r) = ^ is derived by comparing the time needed for correlations 
to span the whole pile and the time needed to create avalanches that span the whole pile (see text). 

This, of course, occurs at a time t x (see eq. [Q, when correlations and avalanches start 
to span the whole 'mountain'. By definition the size of such an avalanche is s x . Since we 
increase the external field at constant rate, t x also is proportional to the number of vortices 
M we must add to the pile to obtain such a 'mountain'-spanning avalanche. However, before 
we create an avalanche of size s x , many smaller avalanches have occurred, adding to the 
number of vortices that we have to add before a 'mountain'-spanning avalanche takes place. 
Thus M is equal to the integral of the size distribution function up to s x , i.e. the shaded 
area in Figure EH This gives 

t x ~ M = sP{s)ds~ s s- T ds ~ s 2 - T (6) 
Jo Jo 
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In the FFS analysis we found that s ~ L , hence we obtain 



- 2-r t D(2-t) 

Combination with eq^ yields the exponent relation 

D(2-t) 



a 

P 



(7) 



(8) 



Equations El and |Hl which were previously derived by Paczuski et a/.[2j, offer an interesting 
possibility: one can calculate the roughness and growth exponents, a and (3, from the 
avalanche properties only. In the table below we compare the values from such an analysis 
with those obtained above from a direct roughness analysis of the surface of the pile: 





a 


P 


from roughness analysis 


0.8(1) 


0.6(1) 


from avalanche analysis 


0.7(1) 


0.5(1) 



Clearly, a good agreement is found, which supports the underlying assumption i.e. that 
SOC theory yields a valid description of the avalanche behavior in this superconductor. 



V. TRANSITION FROM NON-SOC TO SOC 



Naturally, one may wonder whether the avalanche and roughness properties in all super- 
conductors are governed by SOC. To investigate this question, we study these properties as 
a function of quenched disorder. As a sample we use a 500 nm thick Niobium thin film on an 
i?-plane (1102) sapphire substrate, which is a superconductor with T c = 9.25 K, but which 
also absorbs hydrogen from the surrounding atmosphere. The hydrogen enters the Niobium 
lattice interstitially and at low temperatures precipitates in small non-superconducting clus- 
ters, that (like other defects) constitute strong pinning centers. The hydrogen gas absorption 
was always done at room temperature by applying a certain hydrogen gas pressure for 1 hour. 
Typical flux penetration patterns are shown in Fig. 

Avalanches are observed during the slow increase of the externally applied magnetic field, 
in similar manner as described above. The avalanches in the hydrogen-free sample are quite 
compact and all of similar size, while those in the hydrogenated sample (e.g. at 18.1 mbar) 
are more fractal and have a large range of sizes. As a consequence, a finite size scaling 
analysis[13] does work for the hydrogenated sample, but not for the hydrogen free sample. 
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FIG. 7: Magnetic flux distribution in a 500 nm thick Nb thin film on an R-plane sapphire substrate 
after zero field cooling to 4.2 K and subsequently applying a field of 6 mT, for different hydrogen 
loading pressures, as indicated. The black area corresponds to the vortex-free phase. Bright regions 
correspond to the vortex-phase. The arrows indicate a branching magnetic flux protrusion, of which 
the branching is increasing with increasing hydrogen content. 

This, together with a significant change in fractal dimension Q], indicates that (at least 
for these samples) SOC occurs only if a minimal amount of static disorder is present. We 
conjecture that this is a more general phenomenon and that also on granular piles SOC 
behavior depends on the possibility that disorder (randomness) can develop on its surface. 



VI. AVALANCHES IN PATTERNED SUPERCONDUCTORS 

If a superconductor is perforated with a square array of tiny holes (1 — 10 fim in diameter 
and designated as 'anti-dots'), the moving vortices tend to be guided in the direction of 
the lattice vector of these holes jl4]]. In this case, avalanches are still possible, but they are 
guided in similar manner, reducing their fractal dimension, see Fig. |HJ However, still a 
power law distribution of avalanche sizes is observed [15|: the guiding apparently does not 
destroy SOC behavior. 



VII. THERMOMAGNETIC AVALANCHES 



Interestingly, in very similar samples, 500 nm thick Niobium thin films on an A-plane 
(1120) sapphire substrate (and without anti-dots), a completely different behavior is ob- 
served, see Fig. M At low temperatures (T < 5.7 K), upon increasing the field after cooling 
in zero field, initially, the vortex density is rather smooth, but somewhat plume-like, similar 
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FIG. 8: Images of a Pb sample with square array of antidots. Fingerlike dendritic penetration 
occurs (a) at [jlqH = 1.2 mT and T = 4.5 K and at (b) hqH = 1.2 mT and T = 5.5 K, treelike 
dendritic outburst coexisting with smooth flux penetration occurs at (c) hqH = 1.5 mT and T = 
6.0 K, and a smooth profile is found at (d) [iqH = 1.5 mT and T = 6.5 K. In the bottom part of 
the images a saw-tooth-like articact from the magneto-optical garnet is visible. The scale bar in 
each figure corresponds to 0.5 mm. 
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FIG. 9: Magnetic flux distribution in Nb thin film on an A-plane sapphire substrate at 40.0 mT 
after cooling in zero field to the temperatures indicated. The scale-bar indicates the local magnetic 
field. 

to what is shown in Fig. EF- After exceeding a certain distance from the edge, however, 
large avalanches occur, which span immediately nearly half of the width of the sample. At 
low temperature (Fig. EK) these avalanches are much more narrow and less branched than 
at higher temperature (Fig. [U). After each avalanche, the 'smooth' front starts again from 
the sample edge and runs over the previous, still existing, vortex landscape. As soon as it is 
close to the previous front, a new avalanche may occur. The avalanches do not reproduce in 
the same spot from one experiment to the next, indicating that the phenomenon is due to an 
intrinsic instability and not to defects in the samples. Recently Aranson et a/.Q developed 
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a nice theoretical model, which was also used as the basis for numerical simulations, which 
very closely resemble the experimental results (similar experimental results were obtained 
before, see refs 3-9 of ref. while similar models were developed by Shantsev et a/.fl^ 



and Baggio et al. 
equations: 



17|). In this model, the system is described by two coupled differential 



e 



V 2 6-6 + 7 j 2 r 



(9) 



^ { dx 
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dx 



d_ 

dy 



dg_ 
dy 



rH 



where 9 is a reduced temperature, j is the current density and r = r (j, 6) the resistivity 
of the superconductor, g is a local magnetization, defined by j = V x g (s) z. For thin film 
samples, like ours, the second equation is highly non-local, which is described by the kernel 
K. The parameter 7 contains the thermal properties of the film and its contact to the 
substrate, while the parameter r is the ratio between magnetic and thermal diffusion times. 
Small values of r make the system more unstable. For a detailed comparison between model 
and experiment see ref. [lsj ]. 



VIII. CONCLUSION 



We find evidence for SOC behavior of avalanches on the vortex landscape of many dif- 
ferent superconductors, provided that enough (quenched) disorder is present. This SOC 
(avalanche) behavior is preserved in the presence of a square array of strong pins (anti- 
dots). For superconductors close to thermal runaway, only large avalanches occur, all of 
similar size (non-SOC behavior). 
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